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Antiferroquadrupolar Order in the Magnetic Semiconductor TmTe 
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The physical properties of the antiferroquadrupolar state occurring in TmTe below Tq = 
1.8 K have been studied using neutron diffraction in applied magnetic fields. A held-induced 
antiferromagnetic component (fe = (|, |, |)) is observed and, from its magnitude and direction 
for different orientations of H, an 0\ quadrupole order parameter is inferred. Measurements 
below Tn ~ 0.5 K reveal that the magnetic structure is canted, in agreement with theoretical 
predictions for in-plane antiferromagnetism. Complex domain repopulation effects occur when 
the field is increased in the ordered phases, with discontinuities in the superstructure peak 
intensities above 4 T. 
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§1. Introduction 

The number of experimental studies devoted to orbital 
phenomena in magnetic materials has grown rapidly dur- 
ing the last five years. Long-range order involving or- 
bital degrees of freedom has now been reported for a 
broad range of materials, from transition-metal oxides 
(V2O3, LaMn03, etc.) to rare-earth or actinide com- 
pounds. In contrast to 3d elements, lanthanide ions 
are characterized by a strong coupling between the to- 
tal spin and orbital momenta, which rules out pure or- 
bital order. Charge order, an important issue in man- 
ganites, is also not relevant for "normal" rare-earth com- 
pounds because the occupancy of the 4/ shell has a fixed 
integral value. On the other hand, it has long been 
recognized 1, 2 ) that interactions between multipole mo- 
ments of the /-electron wavefunctions — especially the 
quadrupole (QP) moments representing the asphericity 
of the charge distribution — can play a significant role at 
low temperature, and even lead to a phase transition into 
a long-range ordered state. 3 ) In analogy with magnetic 
systems, this state is termed "ferroquadrupolar" when 
the order parameter is uniform, "antiferroquadrupolar" 
(AFQ) when its periodicity differs from that of the lat- 
tice. 4 ) 

A few members of the second class have been known 
for as much as 20 years (CeB 6 , 5 ) PrPb 3 6 )), but more ex- 
amples have been discovered recently (DyB2C2 7 ') as a 
result of an intensive experimental effort using a broad 
spectrum of techniques. The fingerprint of the AFQ 
transition may be elusive in the magnetic susceptibility, 8 ) 
but clear anomalies are observed in the specific heat 5,8 ) 
or the elastic constants. 8 ' 9 -* In most cases, the complex- 
ity arising from the tensor character of the QP moment 
operator can be disentangled only through single-crystal 
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diffraction measurements. 

In this context, x-ray scattering has the unique ad- 
vantage that it can probe the asphericity of the 4/ 
charge distribution without the need to apply a mag- 
netic field. AFQ order was detected recently by this 
technique in several compounds (NdMg, 10 ) CeBg, 11 ' 12 ' 
UPd3, 13 - ) DVB2C2 14 ' 15 **). However, such experiments are 
not straightforward because the measured intensities are 
very weak in the case of conventional Thomson scatter- 
ing, and those obtained in resonant experiments are dif- 
ficult to interpret quantitatively. Neutrons, which have 
no direct interaction with QP moments, can nonetheless 
disclose the "hidden" QP order through its effect on the 
magnetic response to an applied magnetic field. This 
was first demonstrated in the case of CeBg by Effantin 
et al., 16 ^ who observed field- induced antiferromagnetic 
(AFM) peaks below the transition temperature Tq, re- 
flecting the existence of distinct easy axes on the two 
AFQ sublattices. That strategy also proved quite suc- 
cessful with TmTe, as reported below. Furthermore, 
neutron provide an ideal tool to study the magnetic or- 
der which sets in at lower temperature and the interplay 
between the QP and dipole order parameters. 

TmTe is a cubic, NaCl structure, magnetic semicon- 
ductor with an energy gap of about 0.35 eV. 17 > Tm 10ns 
are divalent (4/ 13 ), and the ground-state multiplet 2 F 7 / 2 
is split by the crystal field into one quartet (Tg) and 
two Kramers doublets fXV) and (r 6 ). Neutron scatter- 
ing results 18 ^ indicate that the overall splitting is quite 
small (about 1 meV) and that the ground state is most 
likely Tg, as proposed previously from elastic constant 
and thermal expansion measurements. 19 ' 20 ) AFQ order, 
with a transition temperature of Tq = 1.8 K was dis- 
covered by Matsumura et al. 8,21 ^ using specific heat and 
ultrasonic experiments. The study of QP interactions in 
this system is of particular interest because, unlike in- 
termetallic compounds such as CeB^ or PrPb3, TmTe 
has a very low carrier concentration at temperatures of 
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the order of Tq, and the role of interactions mediated by 
conduction electrons 22 -* can be considered negligible. On 
the other hand, no sizeable lattice distortion indicative of 
magnetoelastic couplings could be detected so far. AFM 
order was reported to occur with a transition tempera- 
ture Tn comprised between 0.23 K (magnetic suscepti- 
bmty 8, 23, 24) and o 43 K ( ncutron diffraction). 25 ) 

In the following, we report a detailed study of the 
AFQ and AFM phases of TmTe by neutron diffraction. 
Section 3.1 presents the observation of the field- induced 
AFM superstructure below Tq, and the analysis of the 
data for different field directions, leading to the identifi- 
cation of the order parameter. These results have been 
described in more detail in previous publications. 26 ~ 29 ) 
Next we set out to discuss the magnetic transition and 
the dipole-ordered phase occurring below T N . Recent 
measurements carried out down to 100 mK are reported 
in Section 3.2. They reveal the existence of a canting of 
the magnetic structure which can be regarded as the con- 
sequence of QP order. The application of large magnetic 
fields at very low temperature produce domain reorienta- 
tion phenomena with dramatic irreversibilities. The final 
section outlines possible directions for future studies on 
this system. 

§2. Experiments 

To observe QP order in TmTe and determine its phase 
diagram, it is necessary to apply large magnetic fields 
and to cool the sample down to T < Tq = 1.8 K. 
Even lower temperatures are required for studying the 
transition into the AFM phase. These conditions were 
achieved by using two different split-coil, vertical field 
cryomagnets with maximum fields of 7.5 and 12 T, re- 
spectively. The former could be equipped with a 3 He- 4 He 
dilution insert yielding a minimum mixing-chamber tem- 
perature of about 100 mK. Thermal coupling of the sam- 
ple to the cold plate carrying the temperature sensors 
was ensured by a gold wire, which proved to be effective 
down to about 150 mK. 

Neutron diffraction experiments were carried out on 
the two-axis lifting-counter diffractometers 6T2 (ther- 
mal beam) and 5C1 (hot source) at the Orphee reactor 
in Saclay. Measurements were made at neutron wave- 
lengths of 0.90 or 2.35 A(6T2), and 0.84 A(5C1) with a 
primary beam collimation of 50'. Second-order contami- 
nations were suppressed by means of erbium or pyrolytic 
graphite filters. For the polarized-neutron experiments, 
5C1 was operated with a Heusler-alloy monochromator, 
providing a polarization Pq = 0.91 at 0.84 A. 

Measurements were taken on a large TmTe single 
crystal (lattice constant a = 6.354 A), prepared by 
induction-melting of high-purity constituents inside a 
vacuum-sealed tungsten crucible. 8 ) From its lattice con- 
stant, as well as from polarized-neutron results reported 
previously, 27 ^ this sample was concluded to be close 
to stoichiometric, with an estimated Tm deficiency of 
0.04 ± 0.01. Different orientations of the crystal with re- 
spect to the magnetic field (H || [110], [001] and [111]) 
were studied in separate experiments. 
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Fig. 1. Intensity of the field- induced magnetic Bragg peak \\\ 
as a function of temperature for different magnetic fields applied 
along [110]. Inset: magnetic form factor associated with the 
induced AFM reflections measured at T = 1.7 K in a field of 
5 T along [111]; the dashed line represents theTm 2 + form factor 
calculated in the dipolar approximation. 

§3. Results 

3.1 Antiferroquadrupolar order 

When TmTe is cooled in an applied magnetic field, 
superstructure peaks corresponding to the wave vector 
k = (typc-II AFM) appear below the critical 

temperature Tq. Intensities measured with if |[ [110] are 
shown in Fig. 1 for different field values. The reflections 
vanish for H = 0, which indicates that the transition 
occurring in zero field is not associated with magnetic 
dipole order. On the other hand, the Q dependence of 
the field-induced intensities approximately follows that 
calculated for the Tm 2+ magnetic form factor (inset in 
Fig. 1). Therefore the observed superstructure results 
from magnetic scattering and cannot be ascribed to a 
distortion of the lattice, which would at any rate require 
unphysically large atomic displacements to account for 
the experimental peak intensities. The observation of a 
staggered magnetic response in a uniform applied field 
is explained naturally by assuming that the transition 
at Tq corresponds to the ordering of 4/ QP moments 
in a long-range structure with the wave vector (5, 5, 5), 
which splits the Tm sites into two sublattices with differ- 
ent local anisotropy axes. Similar results have been ob- 
tained for H || [111]. For H \\ [001] extra peaks are also 
observed but their intensities are much weaker. Phase di- 
agrams drawn from the neutron data for the three field 
directions are in general agreement with those derived 
previously from the specific heat. In particular, the crit- 
ical temperature is strongly enhanced by the external 
field for all symmetry directions, and the QP transition 
line closes up above 5 T for H || [001]. 

To obtain more detailed information on the AFQ order 
we have refined the induced AFM structure for H | [110] 
and H \\ [111]. For both orientations, the AFM Fourier 
component was found to be directed along one of 
the three two-fold cubic axes perpendicular to the wave 
vector k. This gives rise to three possible magnetic do- 
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Fig. 2. Schematic representation of the AFQ structure showing 
the A and B sublatticcs associated with the wave vector fci (left), 
and the field-induced uniform (mp) and staggered (mj^p) com- 
ponents for H || [110] (right). 

mains, denoted "S'-domains" , whose populations can be 
derived from the refinement. From the field dependence 
of the fc-domain populations (see below) the structure 
is expected to be single-fc, and the the AFM moment 
on the Tm sublattices is thus equal to im*,. The re- 
sulting structure for H | [110] is represented in Fig. 2 
for the wave vector (5,5,5)- The AFM component 
{m k = 1.6 [1b f° r H = 5 T) is drawn perpendicular to 
both H and fc because this corresponds to the S'-domain 
found to prevail in high fields. Also shown is the uni- 
form magnetic component along H (m = 2.60(3) Hb 
for H = 4 T from the polarized- neutron results). For 
H || [001], a similar refinement was not possible because 
the experimental intensities were not strong enough. 

Using the symmetry relations between the AFQ order 
parameter and the field-induced dipole moments derived 
previously for CeB 6 , 30 ^ it can be inferred from the above 
results that the primary QP order parameter is most 
likely O 2 . 26 ' A detailed study of multipole interactions 
in TmTe has been carried out by Shiina et a/. 31 ) using 
a classical mean-field model, in which the Tm sites are 
considered to form four uncoupled simple-cubic (sc) sub- 
lattices. Their results indeed confirm that, if the crystal 
field is chosen with an easy axis along [001], the experi- 
mental data can be explained consistently in terms of a 
T 3 (tetragonal) QP order parameter. The situation for 
H || [110] is depicted schematically in Fig. 2. 

The variations of the AFM peak intensities with the 
external field reveal interesting domain repopulation ef- 
fects. In low fields, the four different fc-domains, as- 
sociated with the wave vectors fci = (5,5,5), k 2 = 

5' — 2)' ^ 3 = (2' — 5' 2") an< ^ ^ 4 = ( — \i 5' 5) are 
found to exist. If H is increased along [111], one single 

domain, corresponding to fci \\ H, grows at the expense 
of the other three. This behavior implies that the struc- 
ture is single-fc. A different situation exists for H || [110] 
as two domains, fci = (5, 5, 5) and k 2 = (5, \, — 5), are 
favored in high fields (Fig. 3). As in the previous case, 
the dominant domains turn out to be those whose k vec- 
tor makes the smallest angle (35 degrees) with the exter- 
nal field. The above model, 31 ) completed by introducing 
an additional electrostatic QP-QP coupling between the 
sc sublattices, 32 ^ correctly predicts the single-fc struc- 
ture, as well as the domain repopulation for H || [111], 
but not that for H | [110]. A more realistic treatment 
of these interactions appears to be necessary. 
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Fig. 3. Intensities of the magnetic peaks ±^ ± ^ ± ^ associated 
with the four fc-domains fci . . . ki as a function of the applied 
field H || [110], measured in the AFQ (a, b) and AFM (c) phases. 
Closed (open) symbols denote measurements taken on increasing 
(decreasing) H . Data for each k vector have been normalized to 
the zero-field value at T = 0.1 K (see Section 3.2). 

The field dependences (H \\ [110]) of the induced in- 
tensities of ±5 ± 5 ± 5 reflections corresponding to the 
four fc-domains are plotted in Figs. 3a and 3b for two 
temperatures within the AFQ phase. At T — 1.9 K 
the zero-field state is paramagnetic and the AFM sig- 
nal therefore appears only above a threshold at about 
0.5 T. At T = 0.7 K, on the other hand, the intensity 
increases steeply from H = 0. In both cases the varia- 
tion takes place in two steps, and the intensities from the 
depleted domains k^ and k^ start to decrease only above 
3.5 T, in connection with a pronounced upturn in those 
from the dominant domains. As a result, a sizeable AFM 
component is induced even in the minority domains be- 
fore the fc-domain repopulation is completed. Another 
point to be noted is that the intensities associated with 
fci and fc 2 differ significantly in large fields, especially 
at T = 1.9 K, presumably because a small deviation 
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from the nominal sample orientation with respect to the 
field makes these domains not strictly equivalent. In con- 
trast, their variations remain very similar up to 3 T. A 
double-fc to single-fc transition could produce this type 
of behavior, but it seems incompatible with the domain 
repopulation observed for H | [111]. Another possibility 
is that the fc-domains essentially retain their initial pop- 
ulations up to 3 T, and that the difference between the 
two types of domains reflects their intrinsic staggered re- 
sponse in the applied field. Detailed data collections and 
structure refinements at several values of H should be 
conducted to clarify this point. As could be expected, 
large irreversibilities at increasing and decreasing field 
are observed in connection with the changes in domain 
populations. This effect becomes more pronounced at 
lower temperature. 

3.2 Dipole moment order 

Various experimental results indicate that a second 
phase transition corresponding to the ordering of the 
Tm dipole moments, takes place at Tn <C Tq. 8 ' 23 ~ 25 ) 
The magnetic structure was identified as type-II AFM 
from the observation of neutron Bragg peaks associated 
with the wave vector k = (i, i, ^). 25 ^ At the time that 
result was reported, the existence of AFQ order was not 
suspected and it was therefore natural to assume the 
magnetic structure to be conventional AFII. However, 
the mean-field theory 31 ) suggests that the AFQ order 
should have sizeable effects on the magnetic phase tran- 
sition: it is expected to favor FM correlations, thereby 
causing a substantial reduction of Tn and, in the case 
of "in-plane " AFM order (ordered dipole moments ly- 
ing within the plane of the 0\ quadrupoles) which seems 
appropriate for TmTe, to produce a canting of the AFM 
structure. This canting consists in the rotation of the 
dipole moments on each sublattice by +<p and — cp re- 
spectively with respect to the local anisotropy axis, and 
should be observable experimentally. 

Our measurements for H = confirm the existence 
of AFII magnetic Bragg peaks. In the upper frame of 
Fig. 4, the integrated intensity of the \\\ reflection is 
plotted as a function of T, together with the previous 
results of Lassailly et al. 25 ^ Despite the use of crystals 
from different origins, the agreement between the two 
sets of data is excellent. A tendency to saturation is ob- 
served in our measurements below 0.15 K, probably due 
to insufficient thermalization of the sample exposed to 
the neutron beam. One notes that the magnetic inten- 
sity persists far above the Neel temperature Tn = 0.45 K 
reported in ref. 25. However, the sharp maximum in the 
T dependence of the Bragg peak width on which the lat- 
ter estimate was based could not be reproduced in the 
present experiments. Our data merely indicate a moder- 
ate monotonic broadening of the rocking curves accom- 
panying the intensity drop as the system is heated up to 
0.6 K. The determination of Tn from the neutron inten- 
sities thus remains elusive. We recall that similar diffi- 
culties were also encountered with other techniques, as 
the ac susceptibility 8, 23: 24 - ) shows a divergence at 0.22 K 
but no anomaly around 0.5 K, and the specific-heat peak 
corresponding to the AFM transition exhibits a complex 
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Fig. 4. Integrated intensities of the \\\ (AFM) and 111 (FM) 
peaks measured in zero field as a function of temperature. The 
data of ref. 25 are plotted in the upper frame for comparison. 

structure with a peak at 0.34 K followed by a shoulder 
near 0.5 K. 8 ) These problems may reflect the existence 
of anomalously strong fluctuations, which have been pre- 
dicted to occur because of competing interactions be- 
tween different multipole moments. 

The lower frame in Fig. 4 represents the integrated 
intensity of the weak 111 nuclear peak as a function of 
temperature. The increase observed below 0.45 K clearly 
demonstrates that a superimposed ferromagnetic compo- 
nent develops below Tn in addition to the AFM super- 
structure, with a magnitude of about 0.4 /xb, as com- 
pared to 2.3 hb for the AFM component. The observa- 
tion of this canting component lends additional support 
to the model of ref. 31. Within the limited accuracy of 
the measurement, no significant magnetic signal seems 
to remain above 0.45 K, in contrast to the AFM compo- 
nent. 

Indications as to the dimensions of the magnetic do- 
mains can be derived from an analysis of the width of the 
111 reflection. At T = 0.55 K in zero field, the inten- 
sity is purely nuclear and the width of the rocking curve 
thus reflects the experimental resolution and the mosaic 
of the crystal. For T = 0.1 K, the magnetic intensities 
were obtained by subtracting out this nuclear signal. The 
results for H = and 6.5 T are plotted in Fig. 5, together 
with the nuclear peak. For the sake of comparison, all 
curves have been normalized because their intensities dif- 
fer by orders of magnitude (nucl:FM[ff = 0]:FM[6.5 T] 
Ri 5:1:150). The graph clearly shows that the magnetic 
peak in zero field is substantially broader then the reso- 
lution. This effect can be ascribed to a finite correlation 
length £ of the FM component, estimated to be of the 
order of 50 A, which probably reflects the existence of 
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Fig. 5. Normalized rocking curves through the 111 reflection. 
Crosses and thick lines: nuclear contribution with Gaussian fit; 
open circles and thin line: zero-field FM component with Gaus- 
sian fit; closed circles: high-field FM component. 
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Fig. 6. Intensities of the magnetic peaks ±^ ± i ± i associated 
with the four fc-domains fei . . . ki as a function of the applied 
field H || [110], measured in the AFM phase (expanded plot of 
the low-field part of Fig. 3c). 

magnetic domains with different directions of canting. 
In contrast, the FM component measured at 6.5 T can 
be superimposed almost perfectly on the nuclear peak, 
implying that £ becomes infinite (to the precision of the 
measurement) when the field becomes large enough to 
reorient the domains. This effect takes place at rela- 
tively low fields, as a narrow peak is already recovered 
at if = 2 T. 

The field dependences (if || [110]) of the AFM peak 
intensities for T k, 0.1 K < Tn is displayed in Fig. 3c. In 
contrast to the data for T > Tn plotted in frames a and 
b, finite intensities exist already in zero field because the 
wave vectors of the AFM and AFQ ordered states are the 
same. The measured values are different for k\ and k 2 on 
the one hand, k$ and k^ on the other hand, owing to dif- 
ferent scattering geometries (Lorentz factor and resolu- 
tion). However, it was checked by comparing integrated 
intensities for different sets of equivalent reflections that 
the zero-field populations of the four fc-domains do not 
differ significantly. The curves in Fig. 3c were thus nor- 
malized to the same value for H = 0, and the same 
scaling factor was also applied to the corresponding data 
plotted in the other two frames. 

One can roughly distinguish three different regimes 
as a function of the applied field. Below 1 T (Fig. 6), 



a general decrease in the AFM intensity takes place, 
with initially a weak maximum near 0.15 T for k$ and 
&4. This region is dominated by dipole moment order 
and may involve some reorientation among 5-domains. 
Such a view is consistent with the specific-heat results, 
in which a magnetic-ordering anomaly could be traced 
up to if « 1 T. At higher fields, the intensities for the 
fci and k 2 domains go through a minimum and start 
to increase. The situation resembles that observed at 
T = 0.7 K, i.e. for pure AFQ order. However, for fields 
in excess of 3 T, one observes a succession of jumps, 
with a positive sign for k\ and k 2l a negative sign and a 
smaller magnitude for fc 3 and k^. Similar discontinuities, 
known as "Barkhausen steps " are commonly observed in 
the magnetization curves of ferromagnets. They origi- 
nate from random irreversible domain-wall jumps as H 
is increased. In Section 3.1, it was indeed suggested that 
most of the repopulation of AFQ domains takes place in 
the range H > 3 T. The observation of such a behav- 
ior in an AFM phase is quite interesting and calls for 
further investigations. Surprizingly, no anomalies have 
been detected in the bulk magnetization for this field 
direction. 24 ^ How domain reorientation responsible for 
sizeable jumps in the AFM field-induced component can 
have no fingerprint in the uniform component thus re- 
mains to be explained. 

§4. Conclusions 

The results obtained in the present work provide a 
comprehensive picture of ordering phenomena in TmTe, 
involving both dipole and QP moments. Qualitative 
agreement is found with the main predictions of a mean- 
field model based on a classical treatment of crystal-field, 
quadrupolar and exchange interactions. In particular, 
the canting of the magnetic moments anticipated for in- 
plane AFM order is observed experimentally below Tn- 
On the other hand, a number of questions remain open. 
Probably the most important one concerns the nature 
of QP interactions in TmTe. As noted at the outset, 
the conduction-elcctron-mcdiatcd, RKKY-typc, mecha- 
nism 22 ) proposed for intermetallic compounds is not ap- 
plicable here. Alternative possibilities are a coupling 
through the lattice, as in Jahn- Teller insulators, 33 ^ or 
a superexchangelike interaction between Tm next near- 
est neighbors separated by one Te atom as proposed in 
rcf. 31. Experimentally, it is not simple to distinguish 
between those mechanisms: in particular, it has been 
argued that AFQ order itself should result in sizeable 
lattice distortions below Tq. 34 ) High-resolution x-ray 
diffraction experiments should be performed to clarify 
this point. Resonant experiments could also allow the 
AFQ order to be probed directly in zero field. This re- 
quires, however, to reach temperatures of less than 1.8 K, 
and to obtain crystals with very small mosaic spread. 

Our present interpretation of the AFQ phase as due 
to an 2 QP order parameter seems consistent with the 
bulk of experimental observations, but some points re- 
quire further examination. In particular, the results for 
H || [001] are not well understood, mainly because of the 
weakness of the signals. Even so, the fact that a finite 
induced AFM component exists seems in contradiction 
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with the symmetry arguments of Shiina et aZ. 30 ) which 
predict that component to vanish. This may indicate 
that the latter analysis, initially performed for CeBg (sc) 
is not readily applicable to fee TmTe because the small 
groups of the k vector (5,5,5) for the two types of Bra- 
vais lattices are different. 3 5 > It was also predicted that 
for this field direction, the external field should stabilize 
the component OS, with respect to 0\, possibly causing a 
transition between two AFQ phases. 31 - ) Experimental ev- 
idence for such a transition is still lacking but it has been 
suggested that it might be responsible for the steplikc 
jump observed in the bulk magnetization for H || [001] 
at very low temperature. 24 ' As to the properties of the 
AFM state, the results reported in this paper reveal quite 
interesting behaviors, especially in strong applied mag- 
netic fields. Additional experiments are required to work 
out a reliable refinement of the magnetic structure in the 
presence of complex domain reorientation effects. 
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